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Core collapse of massive stars resulting in a relativistic fireball jet which breaks through the stellar 
envelope is a widely discussed scenario for 7-ray burst production. For very extended or slow rotating 
stars, the jet may be unable to break through the envelope. Both penetrating and choked jets will 
produce, by photo-meson interactions of accelerated protons, a burst of > 5 TeV neutrinos while 
propagating in the envelope. The predicted flux, from both penetrating and choked jets, should be 
easily detectable by planned 1 km^ neutrino telescopes. 
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The leading model for gamma-ray bursts (GRB) in- 
volves a relativistic fireball, where the 7-rays are pro- 
duced by synchrotron or inverse Compton radiation from 
Fermi accelerated electrons in optically thin shocks. The 
ultimate energy source of the fireball is thought to be the 
gravitational energy release associated with temporary 
mass accretion onto a black hole, resulting either from 
a compact merger |^ or from the collapse of a massive 
star Q (see [|| for a review). All the GRB so far local- 
ized with the Beppo SAX satellite [Q belong to the class 
of long bursts ||] with durations ^ 2 s, and there 
is increasing evidence that these occur in star-forming 
galaxies For this class at least, comprising ~ 2/3 of 
all bursts, it is widely assumed that a massive, collapsing 
star ( "coUapsar" ) is the progenitor. The shocks produc- 
ing the 7-rays must occur after the fireball has emerged 
from the stellar envelope. The preferred escape route is 
along the centrifugally lightened rotation axis, and the 
stellar pressure tends to collimatc the fireball into a jet. 
The shocks should also lead to Fermi accelerated rela- 
tivistic protons 0,^, and to neutrinos with e^, > 10^'' eV 
via interactions with the ~ MeV 7-rays [^,|l^ (see ||l^] 
for a review). 

In this Letter we show that while a jet is making its 
way through the star, its rate of advance is slowed down 
in a termination shock to a Lorentz factor much lower 
than it has at injection or after emergence; and additional 
internal shocks are expected in the jet interior, inward of 
the termination shock. The latter can accelerate protons 
to ^ 10^ GeV, which interact with X-ray photons in the 
jet cavity leading to electron and muon neutrinos (and 
anti neutrinos) with energies £1/ ^ 5 TeV. These neutrinos 
appear as a precursor signal, lasting for time scales ^ tens 
of seconds prior to the observation of any 7-rays produced 
outside the star by a coUapsar-induced GRB. 

Furthermore, in a significant fraction of the massive 
stellar collapses, the jet may be unable to punch through 
the stellar envelope However, the TeV neutrino sig- 
nal from such choked jets should be similar to that from 



"successful" jets which do break through and lead to ob- 
servable GRB. This may provide a means of detecting 
and counting such choked-off, 7-ray dark collapses. 

The TeV fluence from an individual collapse at cosmo- 
logical distance z ^ 1 implies 10~^ — 10 upward moving 
muons per collapse/burst in a 1 km^ detector. The sig- 
nal from an individual collapse with associated electro- 
magnetic signals (and possibly from rare, energetic dark 
collapses), are well above the atmospheric neutrino back- 
ground, and the total integrated flux from dark collapses 
should also be detectable. The AMANDA experiment 
Jl3| , p^ may soon be able to provide relevant limits on 
the total integrated rate of such dark collapses. 

Jet propagation inside a massive progenitor envelope. 
Numerical simulations of coUapsar models leading to 
black hole formation have been performed for a range 
of stellar progenitors and initial conditions. The charac- 
teristic progenitor structure is that of an evolved mas- 
sive star, with ^ "^Mq Fe core of radius ^ 10^ cm and 
~ 8Mq He core extending to ~ 10^^ cm, in some cases 
surrounded by an H envelope extending to 10^'^ cm (a 
red supergiant), while in others the envelope has been 
(largely) lost. While the details remain uncertain, pre- 
liminary calculations suggest that a relativistic jet can be 
launched along the progenitor rotation axis [|l2|,^ , pow- 
ered either by thermal annihilation or MHD stresses 
coupling to the black hole and to the debris disk falling 
back onto the hole. The jet life time is limited by the gas 
fall-back time onto the black hole, and should be compa- 
rable to GRB durations, ^ 10^ s. (Similar pair or MHD 
outflows may also arise from collapses leading to highly 
magnetized pulsars ]l^ , [l7t ). 

When 7-rays are observed, the required isotropic 
equivalent luminosities are Liso. lO^^erg/s and the in- 
ferred jet Lorentz factors are Tj ^ 10^. This is the fi- 
nal value outside the star, which may be representative 
of the intrinsic injection Lorentz factor. As it advances 
through the star, the jet drives a bow shock ahead of 
it. The jet is capped by a termination shock, and a re- 
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verse shock moves back into the jet, at which the rela- 
tivistic jet plasma is decelerated. Both the shocked jet 
plasma and the shocked stellar plasma advance with a 
jet head Lorentz factor Th <C Tj [Q. In the rest frame 
of the shocked, decelerated plasma, the fast jet moves 
with a Lorentz factor « Tj/2Th- During propagation 
in the Fe and He cores, the high density ahead of the jet 
implies deceleration of jet plasma to subrelativistic veloc- 
ity. However, beyond the He core edge at r ~ 10^^ cm, 
the density in an H envelope drops to p ^ 10~^gcm~'^, 
and the jet can accelerate to relativistic velocity (c.f. 
also Balancing the energy density behind the for- 

ward shock, 4r^pc^, with that behind the reverse shock, 
A(T j /2T fi)'^njmpc'^ , where the jet proper proton number 
density nj = Liso/(47rr^r^mpC'^), we find 
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where Liso = 10^^L52erg/s, r — 10^^ri2 cm is the shock 
radius, and p = 10~^p_7g cm^^. 

The jet head propagates at mildly relativistic veloc- 
ity Th up to a typical super-giant H envelope radius 
r ~ 10^'^ cm (and if it makes it through entirely, in 
the even lower density beyond it would increase back to 
- lO^Fa). Envelopes larger than i?, > 10"i?i3 cm, or 
denser cores, e.g. due to slower rotation rate, could how- 
ever lead to a stalling of the jet before it breaks through 
the star, if the jet crossing time is longer than the jet 
lifetime (as argued by, e.g. [T^Jl^ ). 

The proper number density n of protons and elec- 
trons in the shocked, decelerated jet plasma is related 
to the proper density nj in the fast moving jet through 
n — 4:(rj/2Tii)nj, and the energy density in the shocked 
plasma is U = 4:(rj/2Th)'^njmpC^. The hot, shocked 
jet plasma drives the bow shock wave that propagates 
forward (radially) and sideways (tangentially) into the 
envelope. For a jet of small opening angle 6 < l/Th 
Ip^ , the proper thickness of the shocked jet plasma shell 
is A ~ Or. This can be seen by balancing the parti- 
cle flux across the reverse shock, « TT{6r)'^c{rj/2rh)nj, 
with the tangential flux of particles leaving the cylin- 
der of height A and radius 9r of shocked plasma, w 
2TT9rA{c/ ^/3)4:{T j /2T h)nj (taking the tangential flow ve- 
locity to be comparable to the sound speed in the shocked 
plasma). Adopting A = 0.29r and 9 = 10"i6'_i, the 
Thomson optical depth of the shocked plasma shell is 
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where F, 



lO^F, 



The reverse shock may start out being collisionless, 
heating the electrons and protons to ~ 10^F2 GeV. How- 
ever, the electrons would lose all their energy on very 
short time scale, by synchrotron and inverse-Compton 



emission, converting a large fraction of the shocked 
plasma internal energy U to radiation. Due to the high 
optical depth, the radiation will thermalize, with an ap- 
proximate black body radiation temperature 
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The corresponding proper photon number density is 
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Thus, the reverse shock is likely to become radiation 
dominated, i.e. the dissipation of jet kinetic energy 
may be mediated by inverse-Compton scattering of the 
electrons, whose mean-free-path is ~ ~ 1 cm, 

where ctt is the Thomson cross section, and electrons 
lose all their momentum in a small number of scatter- 
ings (scattering of a thermal photon boosts its energy to 
^ 10(F2/F/i)^ MeV in the shocked plasma frame). 

The mean-free-path of thermal photons propagating 
back into the jet, in the decelerated plasma frame, is 
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Thus, the thermal photons do not penetrate much into 
the jet. If the reverse shock is indeed radiation domi- 
nated, the shock thickness would also be of order l^. 

Photo-meson Interactions and Neutrinos. In addition 
to the bow and reverse shocks moving with Th <^ Tj 
in the star, internal shocks in the pre-deceleration jet 
moving with Tj are expected to occur at radii 



T^jcSt = 3 X 10"F^(5t_3 cm. 
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which are smaller than the termination and reverse shock 
radius. Here 6t is the variability timescale of the injected 
relativistic outflow, whose minimum value s ~ lO^^St-3 
s mjl^. While the details could depend on the ratio of 
MHD to pair and baryon energy, it is commonly assumed 
that internal shocks are mildly relativistic in the jet 
frame, and protons will be Fermi accelerated to a power 
law distribution which approximates dup/dep oc e~ ^ (see 
jm for detailed discussion). 

As they approach the reverse shock, high energy pro- 
tons may produce pions in photo-meson interactions with 
thermal photons of energy ~ ~ 4 keV, provided 
their energy (measured in the decelerated plasma frame) 
satisfies epe^ > 0.3GeV^ i.e. for ep > 10^ GeV. The 
maximum energy to which protons may be accelerated 
in the internal shocks is limited by synchrotron losses. 
The acceleration time, in the jet frame, to observed en- 
ergy e'p = F/iEp, corresponding to energy « in the 
decelerated plasma frame, is approximately 1^, 11 



t — = 

c TjcB 



= 10 



-6.5 



TlThSt^ 



(^52^5,-2)1/ 



YT^Ep.iss, (7) 



2 



while the (jet frame) time scale for synchrotron losses, 
which dominate at high energy, is 
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Here, Cp = lO^^Cp^s eV, Rl is the Larmor radius and 
= 10~^C_B,-2 is the magnetic field equipartition frac- 
tion, Anr'^cr'jB'^ /Stt = ^siiso.- Thus, protons may be 
accelerated to energies well above the pion threshold. 

Protons at the threshold energy, for which the ther- 
mal photons are at the A-resonance, lose w 20% of their 
energy in each interaction. Inside the jet, the thermal 
photon density drops exponentially with the distance x 
from the reverse shock, nj{x) cx exp{—x/lj), and reso- 
nant protons will lose all their energy to pion production 
at a distance x from the termination shock where the 
photon density satisfies 
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Here, = 5 x 10~^*cm^ is the cross section at reso- 
nance. Typically, resonant protons and neutrons (which 
are produced with tt+'s in proton photo-pion interac- 
tions) will lose their energy at x/l^ « 10, while higher 
energy nucleons, for which lower, will lose their 

energy to pions at shorter distance. 

The pions produced in photo-meson interactions lose 
energy by inverse-Compton scattering on thermal pho- 
tons. The lowest energy pions produced are char- 
acterized by Lorentz factors comparable to those of 
the resonant protons, w 10^(Tr/4keV)~^ in the de- 
celerated plasma frame. Inverse-Compton scatter- 
ing of ~ 4 keV photons by pions of higher energy 
is in the Klein-Nishina regime. Thus, the inverse- 
Compton energy loss time of high energy pions is Tie ~ 
[3n!;crT(me/m^)2c/8]-i(e^/10i3eV)(Tj4keV). The ra- 



tio of Tic to the pion life time is 
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We therefore expect « 1/4 of the energy lost by protons 
to pion production to be converted to muon neutrinos 
(w 1/2 of the energy is lost to neutral pion production, 
and ~ 1/2 of the charged pion energy is converted by 
decay into muon neutrinos). 

If the termination shock is coUisionless, rather than ra- 
diative, protons may be accelerated to energies exceeding 
the photo-meson threshold also in the termination shock. 
However, it is clear form Eq. ( p^ ) that inverse-Compton 
losses of the pions in the high density photon field pre- 
vent the production of very high energy neutrinos from 
protons accelerated in the termination shock. 

The Vfji,{v^, Ve, Vf.) from pion decays due to protons ac- 
celerated in internal shocks will have a typical energy 



5 — 10%ep. We therefore expect a flat power per decade 
neutrino spectrum H, e^dn^/dei, oc e°, at 
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We note, that since photons propagating backward into 
the jet are emitted from a region of shocked, decelerated 
plasma with small optical depth, the spectrum of photons 
with which protons interact may differ from thermal. In 
particular, a non-thermal extension of the photon spec- 
trum to energies 3> may lead to an extension of the 
neutrino spectrum to energies < 1 TeV. A detailed cal- 
culation of the photon spectrum is, however, beyond the 
scope of this Letter. 

The jet isotropic-equivalent luminosity, iiso. ~ 
lO^^erg/s, and energy, i^iso. ~ lO^'^ erg, are inferred 
from 7-ray observations, and reflect the energy of elec- 
trons accelerated to high energy by the internal, mildly- 
relativistic coUisionless shocks. Protons are generally ex- 
pected to be accelerated in such shocks with higher ef- 
flciency (see e.g. for review). We therefore assume 
that the accelerated proton energy is at least compara- 
ble to that inferred from observations to be carried by 
accelerated electrons (see also |^-^,^). A large fraction 
of this accelerated proton energy will be converted, while 
the jet is inside the envelope, to pions. The envelope is 
transparent to TeV neutrinos, and the expected TeV neu- 
trino fiuence at Earth from a jet at a luminosity distance 
Dl ^ lO^s cm (redshift z - 1) is 
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where Ei^o. = lO^^Ess erg. 

The probability that a muon neutrino will produce 
a high-energy upward moving muon in a terrestrial de- 
tector is approximately given by the ratio of the muon 



range to the neutrino mean free path [gO|, P^^, w 1.3 x 
10~^€.TcV' with /3 = 2 for e^^^xoV < 1 and /3 = 1 for 
(10) Ei/.TcV > 1- For a flat neutrino spectrum, this implies that 
the number flux of muon induced neutrinos can be re- 



lated to the neutrino energy flux F^, by Jp — [Pq,/ eufl)Fu, 
with Po/e,.,o = 1-3 X 10-^TeV"^ Using Eq. (|l|), the 
average number of upward muon detections for an indi- 
vidual burst is 
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Implications. The typical angular resolution of 
planned neutrino telescopes at TeV energies (see |l^] 
for a review) should he 9 ^ 1 degree. The 
atmospheric neutrino background flux is '^u,bkg 



10- 



-2.5 



cm ^s ^ster ^ , implying a number of > 1 TeV 
background events N^kg ~ lO~^(0/deg)^i2km~^ per an- 
gular resolution element over a burst duration of 10^<2 s. 
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Thus, the neutrino signals from individual collapses at 
z ~ 1 would stand out with high statistical significance 
above the background. 

The rate of 7-ray-detectable GRB events, resulting 
from successful breakthroughs of jets beamed towards 
Earth, is ~ lO'^/year j2l]. The predicted neutrino signal, 
~ 0.2/km2 muon events/burst at e^^ > 5 TeV [Eq. @], 
would precede by tcross ^ tj = 10^i2 s the GRB trigger. 
Its detection would be a test of the coUapsar hypothesis 
for long GRB, and the duration of the neutrino precur- 
sor would constrain the dimensions of, and jet advance 
speed in, the stellar envelope. In this case the absence (or 
a much reduced duration) of i^-precursors in short GRB 
would be a test of the view that long and short GRB 
arise from qualitatively different progenitors. 

The signals from both emerging and choked-off jets 
in equation are calculated for "mean" events, i.e. 
for E,so ^ lO'^^^^ss erg, - lO^Es and z - 1. The 
total fireball jet energy and bulk Lorentz factors will, 
however, be distributed around these mean values, and 
fluctuations around the mean in these quantities and in 
the distance will lead to detections dominated by rare, 
energetic nearby events, as argued for GRB Thus, 
one could expect a significant number of neutrino bursts 
(few/year) with A^^ > 10/km^ muon events/burst. These 
should be easily distinguishable above the atmospheric 
neutrino background at z ~ 1, whether coincident with 
observed GRB or not. 

The number of coUapsars with Earth-pointing choked 
jets, while model-dependent, may be in excess of the 
number of successful collapsar GRB (where jet punch- 
through occurs and 7-rays are detected), ^ 10'^ /year. If 
jet punch-through requires fast core rotation rate, e.g. 
from a massive star merger with a binary companion in 
a fraction of cases, choked coUapsars could greatly ex- 
ceed successful ones. The number of massive collapses 
may, e.g., be comparable to that of type II supernovae 
(SNe II), ~ 3 X 10^^/year/galaxy. For a galaxy den- 
sity no ^ lO^^/Mpc"^ within a Hubble radius ~ 3 Gpc 
{z ~ 1) the total rate of SNe II is Usn 1 s"^ If all 
SNe II lead to jets beamed into a solid angle ^Ij ^ 10~^47r 
(e.g. |24yi(i|j25|] ), the effective rate beamed towards Earth 
would be ^ lO'Yyear. From equation (12), the vv energy 
flux per logarithmic energy interval from all directions 
would be ti^u ~ 10~^GeV/cm^s sr, where we divided 
by \n{e^^,nax I ^y,min) ~ 10. This signal is not far below 
the experimental upper bound on the diffuse neutrino 
flux recently established by the AMANDA experiment, 
e^'&i^ < lO^^GeV/cm^ssr iQ. An increase in neutrino 
telescope exposure will therefore allow in the near future 
to put relevant constraints on the frequency of 7-ray dark 
i^-bursts from choked jets, and hence on the frequency of 
core collapses at the poorly known high end of the stellar 
mass distribution. (Unfavorably aligned emerging jets 
would also appear 7-ray dark, but would be neutrino- 
dark as well). 



The TeV neutrino signals preceding the 7-rays in GRB, 
as well as those from choked, 7-ray dark jets, would differ 
significantly from other neutrino signals expected during 
and after the 7-ray phase of GRB. In the latter, one 
expects ^ 100 TeV neutrinos |^ from internal shocks 
well beyond the stellar envelope accelerating protons that 
interact with MeV photons, or ^ 10^^ eV neutrinos p^ ] 
from external shock protons interacting with UV photons 
even further out, and one also expects 2-5 GeV neutrinos 
p7|^-p9[ from inelastic nuclear collisions when neutrons 
decouple from protons. For the flat e^'^v neutrino spectra 
considered here, the number of precursor TeV neutrino 
events per burst is also larger by at least one order of 
magnitude, relative to those at ^ 100 TeV and at GeV 
expected during and after the 7-bright phase of GRB. 
Such GRB precursor or 7-ray dark TeV neutrino burst 
signals may therefore be the likeliest targets for early 
detection with planned km^ experiments . 

If, as suspected, the first generations of stars formed in 
the Universe are much more massive than those forming 
now, then at redshifts z 5 jet break-through collapses 
may be rare, and 7-ray dark choked jet collapses leading 
to TeV neutrinos may be much more common than today. 
Without an emerging jet, the X-ray, optical, IR and radio 
afterglows typical of GRB would also be absent. Such 
neutrino bursts, which at z 5 produce a ^ 1 TeV 
fluence smaller by only a factor of « 10 compared to z ~ 
1 bursts, would therefore be a prime tool for investigating 
both successful and failed GRB, as well as massive star 
formation at such high redshifts. 
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